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Abstract 

The local-density approximation (LDA), together with the half-occupation (transition state) is 
notoriously successful in the calculation of atomic ionization potentials. When it comes to extended 
systems, such as a semiconductor infinite system, it has been very difficult to find a way to half- 
ionize because the hole tends to be infinitely extended (a Bloch wave). The answer to this problem 
lies in the LDA formalism itself. One proves that the half-occupation is equivalent to introducing 
the hole self-energy (electrostatic and exchange-correlation) into the Schroedinger equation. The 
argument then becomes simple: the eigenvalue minus the self-energy has to be minimized because 
the atom has a minimal energy. Then one simply proves that the hole is localized, not infinitely 
extended, because it must have maximal self-energy. Then one also arrives at an equation similar 
to the SIC equation, but corrected for the removal of just 1/2 electron. Applied to the calculation 
of band gaps and effective masses, we use the self-energy calculated in atoms and attain a precision 
similar to that of GW, but with the great advantage that it requires no more computational effort 
than standard LDA. 

PACS numbers: 71.15.-m 31.15.-p 71.20.Mq 
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I. INTRODUCTION 



The well known Density- Functional Theory (DFT)i is an approach to the theory of elec- 
tronic structure in which the electron density distribution, rather than the many electron 
wavefunction, plays a central role. The practical applications of DFT are based on approx- 
imations for the so-called exchange-correlation potential which describes the effects of the 
Pauli principle on the many-electron system. If we had the exact exchange- correlation po- 
tential we could solve the many-body problem exactly for the ground state. Although the 
potential is unknown, approximations are made. The most common is the so-called Local- 
Density Approximation (LDA) which locally uses the exchange-correlation energy density of 
a homogeneous system. The LDA to the Kohn and Sham DFT- is still one of the most reli- 
able methods for condensed matter calculations, having successfully predicted and explained 
a wide range of ground state properties in solid state physics and chemistry^. Lately, but 
very slowly, it is being progressively abandoned in favor of the many GGA (Generalized Gra- 
dient Approximations)^. However, while LDA and GGA have predicted many ground state 
properties with good accuracy, the electronic properties such as band gaps are significantly 
smaller than experiment. These discrepancies are caused by the lack of the discontinuity 
of the exchange-correlation potential^ in going from the valence to the conduction band. 
Several methods for overcoming these limitations have been proposed. One of them is 
the GW approximation, in which one considers the energies of quasiparticles and calculate 
the electron self-energy in terms of perturbation theory.- 1 ^ This procedure has been quite 
successful, achieving good accuracy, but it goes beyond the DFT. Other procedures were 
also proposed, among them we can mainly cite the Self Interaction Correction (SIC)£, the 
atomic SIC applied to solids (ASIC)^^, which is perhaps the procedure closest to ours, hy- 
brid functional^, screened exchange (SX-LDA)ii, the so-called exact-exchange approach^, 
the well known LDA+U— , the work of Liberman^, and others. Most of these approaches 
are computationally very demanding, which prohibits their application to large systems of 
atoms. 

Slater half-occupation scheme^ 1 ^ 1 ^ was very successful for valence states. One example 
is that one could obtain energies which were comparable to the experimental ionization 
energies^, though, at that time, good spin-polarized exchange- correlation approximations, 
as those based on Ceperley and Alder, did not exist.- In order to illustrate to the reader 



2 



the quality of the results that can be obtained, we present in Table d the first and second 
ionization potentials of 12 atoms, measured and calculated with LDA with 1/2 occupation. 

TABLE I: First and second Ionization Potentials (IP) for some atoms (eV). These results were 
obtained with spin-polarization but assuming spherical charge densities for ions and atoms. We 
used a code originally written by S. Froyen, modified by N. Troullier and J. L. Martins, and 
modified and maintained by A. Garcia. 
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Though the precision of the calculated results shown in the Table is much better than the 
precision one reaches in the calculation of band gaps, either by LDA or GGA, it has been 
difficult to find a way to make the ionization of 1/2 electron in extended systems as crystals. 
Of course the problem is that a crystal is described by means of Bloch waves and removing 
the population of just one Bloch state is of no consequence. In this paper we present a 
solution to this problem. We are specially concerned with the calculation of band gaps in 
semiconductors, for which we obtain calculated results that compare very favorably with 
experiment, and are not computationally demanding. We report the results for fourteen 
semiconductors, including the group II- VI, III-V and IV. Our method is inspired in the 
LDA and in the half ionization but, at some point, it has to be postulated. The quality of 
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the results and the ease with which they are obtained show that our assumptions are very 
good. Now we develop our method which could be properly named as LDA-1/2. 



II. LDA AND HALF-IONIZATION IN SOLIDS 

Accepting the LDA as a valid approximation to the density functional (DF), the Janak 
theorem^ follows 

8E 

Wa =eMa) (1) 
where E is the total energy of the system, function of the occupation f a of the one-particle 
Kohn and Sham state a. It is a well known fact that the eigenvalue e a (f a ) is almost precisely 
linear with the occupation f a .— Then integrating 



dfa 

between the ground state (f a = 0) and the ion (f a = —1) obtain 

E(0) - E(-l) = e a (-l/2) = - Ionization Potential. (2) 
Taking another derivative 



where 



S. = / / ^-f^n + M , M f )l -^ r) n a( r <) (4) 



\r-r'\y df a 2JJ 5n(f)5n(r')y df a 

is named "self-energy", because of the first term in the right. In LDA the functional deriva- 
tives become common derivatives times delta-functions. We maintain the functional deriva- 
tive notation because the final formulae can have extended use. Because of the linearity of 
^a(fa) —we may write 

e Q (-l/2) = e Q (0) - S a (5) 

and 

E(0) = E(-l) + e a (0)-S a (6) 

Eq. [6] is quite surprising in its simplicity. The Eq. is telling us that, to restore the ground 
state, with total energy E(0), from an ion with a hole at state a we add an electron whose 



energy is the eigenvalue e a (0) minus the hole self-energy. The self-energy is large when 
the function is much localized as an atomic wavefunction, and is small and zero when it is 
much spread as a Bloch function. Since the energy of the restored ground state must be a 
minimum the hole self-energy must be a maximum. Thus the hole should be representable 
by a very localized wavefunction. So far, this is the first time we see a demonstration of the 
hole localization, though this proof is based on an approximation (LDA) to the DF theory 
and on the linearity assumption. Of course we cannot say that the localized hole state is 
truly stationary, specially if its energy is inside the band continuum of the Bloch states, into 
which the localized hole would be scattered. 

The self-energy may be thought as the Quantum Mechanical average of a "self-energy 
potential" Vs(r) such that 

S a = J d 3 rn a (r)V s (7) (7) 

where n a = ^ip a , 

V s (f) = I d?r ' + \j£r •-J^- FJ y a (f ') ( 8) 



3 ^pfp^- I r 3 , 5 2 E XC v dnp{r ') 



8n{r)8n{f ') p df a 

that depends on the state a. From now on, in Eqs. such as[8l we will not write the last two 
terms, those depending on the derivative of the wavefunctions with respect to the occupation 

fa- 

To derive Eq. [6] we assumed linearity, aside from the Janak theorem. The linearity results 
when the Kohn and Sham eigenfunctions of the ground state are equal to those of the ion, 
which is correct to a large extent.— Coherently we may neglect the last two terms in Eqs. 
such as [H and make the difference E(0) — E(—l), of two minima, an extremum. Then we 
minimize (extremize) e a (0) — S a as suggested by EqEJ To do so we must write a variational 
expression that, upon minimization (extremization), leads to a differential Eq. for the hole 
wavefunction ?p a (r). At this point we must set clearly what we imply by the term "hole". 
From Eq. we see that we are adding an electron to a hole state of the ion. The hole state 
might be in a valence or conduction band, the only requirement being that the state is empty 
in the ion. Thus by a "hole" we mean "an electron filling an empty state" or u a particle 
excitation" . This particle excitation may be in the valence or in the conduction band. 
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The variational expression must be such that upon calculating it with the solution ip a {r) 
of the differential Eq. it returns back the average e Q (0) — S a . Then we write it as 



j \r — n\ J \r — r \ ono{r) 



4><* (9) 



where n (f) does not include the hole wavefunction ip a and it is the number-density of the 
ground state. In this case, the wavefunction tp a is to be interpreted as the hole state to be 
filled with an electron to restore the ground state. 

Performing the extremization and then inserting EqlHlwe obtain the equation below with 
the top entries of the brackets {} 

-V 2 " 2 £ A + 2 / ^Qj-A ' + Pf- (10) 
j \r — rj\ J \r — r \ ono{r) 



{ } / p — tfV -{ } —-———-ip a {r ) ip a (r )d 6 r 

2 J \r — r \ i J dno{r)dno{r ') 



If we insert Eq. |8] before extremization we obtain the Eq. [10] with the entries at the bottom 
of {}. In this latter case, the next to last term of the operator is exactly the term in 
the SIC equation^. Its effect is to exclude the electron being considered from the Hartree 
interaction. The last term in the operator, the exchange- correlation term, is very different 
from the corresponding SIC term, since it depends on the whole density of the system and 
not only on the density of the a state. The SIC equation, Eq. [10] with the bottom entries 
in {} is not what we want because the eigenvalue X a = e a — 2S a and not X a = e a — S a as 
the half- ionization requires. It is worth mentioning that, in the calculation of band gaps, 
SIC over corrects and halving it seems to be a better procedure.- 

Except for atoms, solving Eq. [TO] is very difficult, either with the top or bottom entries 
in {}. One important problem is that the solutions of Eq. [10] are not orthogonal. The SIC 
solution for atoms is used in the ASIC method, which is excellently reviewed in ref. 0. In 
our case we proceed differently. We introduce a parametrized self-energy potential and use 
a variational expression that is an extremum for variations in the parameter (s). The first 
question to be answered is whether it is possible to define a unique self-energy potential that 
is state- in dependent. In Table [IT] we show a study of how the atomic self-energy of many 
states vary with the assumed self-energy potential. One sees that for s and p orbitals, the 
self-energy does not vary much whether it is calculated with s, p, and even d self-energy 
potentials. Of course the self-energy potential that we will use is the one corresponding to 
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the atomic orbital dominating the crystal energy bands around the gap. In extreme cases 
we can define a self-energy potential that is angular momentum dependent. 

I 

V s = Y, V sAr) E \l,m><l,m\ (11) 

l m=—l 

This possibility was explored in the case of diamond, as shown in the discussion of our 
results. 

TABLE II: Atomic self-energies (eV), Eq. for some third-row atoms and different valence elec- 
trons, calculated with self-energy potentials (Eq. [TBj) derived from half- ionization of different states. 
The fact that the self-energy of a given state has no important dependence on the self-energy po- 
tential, that is the entries at each column for a given atom do not differ much, shows that, for 
practical purposes, one can neglect the state-dependence of the self-energy potential. In the table 
below, only the self-energies of d-states seem to be much dependent on the way the self-energy 
potential is derived. 
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Thus, assuming a self-energy potential that is state-independent, we use the following 
variational expression.— 1 ^ 

E[n, v, p] = K[n] - J V[p]p + \ J V[p}p - J V s p 
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+ J v(n-p)+E xc [p] (12) 

where p is the proton number- density, n = J2p fp^pi^p is the electron number-density made 
out of the squares of the wavefunctions, fg being the occupation numbers, v is the Kohn- 
Sham potential, p is the model number-density, Vs is the given parametrized self-energy 
potential, K and E xc have their usual meaning of kinetic and exchange-correlation of the 
Kohn-Sham DFT. The functional V[p] is defined as 

VW)\ = 2 / p^A' (13) 
J r — r 

The functional E is an extremum for variations in any of the three functions n, v, and 

P ■ 

1. 5E/5v = leads to p = n 

2. 5E/5p = leads to 

v = -V[p] + V[p] -V s + 5E xc /5p (14) 

3. 5E/5n = leads to Schroedinger Equations with potential v and eigenvalues e a , 
rewritten as. 



j |r — r/| «/ \t — t | op(r) 



ip a (r) = e a ip a {r) (15) 



Using Eq. dH to determine p for given v and Vs and solving the Schroedinger equations, 
find 

E = EfM-lf y\p)P + EM -jp^f (16) 

It must be understood that both n and p are number densities of iV electrons, not N — 1/2, 
but the eigenvalues correspond to a situation where 1/2 electron is removed, if Vs is well 
chosen. 

We want to find band gaps by taking the difference of total energies due to different 
occupations f a . Maintaining the Kohn and Sham potential v and the model number- density 
p, solution of Eq. Q3J the band gap becomes a difference between eigenvalues e a . Now, 
because the total energy is a variational functional, that is an an extremum for variations 
in v, resulting from variations of the self-energy potential Vs, one should look for extreme 
eigenvalue differences. 

= (17) 
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III. THE LDA-1/2 METHOD 



Consider the case of an atom. We first prove that the self-energy potential is given by 

V s ~-V(-l/2,r) + V{0,r) (18) 

namely the difference between the all-electron potentials of the atom and of the half-ion. 
We begin by writing the potential difference as 



V(-l/2,r) - V(0,r)= / ""V » 1-2* +2 />r< "<«"> ■ 



xc 



o <9/ a [ r J |r — r'| 5n(r) J 

o 1/2 *4 2 / d3r '^ + / d3r '^SH MF,) } (19) 



Jo J \r — r'\ J on(r)bn{r ') p of a 

or, for a certain value of f a in [—1/2, 0] 

-V(-l/2,r) + \/(0,r)= fd 3 r' " a ^\ +- I d z r 1 —^^—n a {f ') 



2 j 5n(r)8n{r ') 

Fig. [T] depicts r-times the self-energy potential for the nitrogen atom, a typical case, for 
degrees of ionization I ranging from 0.5 to —0.2. Observe that the ratio 

V(fa,r)-V(0,r) 

fa 

has a very poor dependence on f a meaning that in Eq. [20] we can take the xc functionals at 
the full occupation f a = 0. Then, comparing Eqs. [HI and [201 our proof is completed. 

We will leave to another paper a discussion on the lone hole solution we can get out of 
Eq. [TUJ Here we are interested in calculating band gaps of semiconductors. For that purpose 
we will repeat the atomic self-energy potential (Eq. I2"U|) in the whole lattice and calculate 
eigenvalues for "hole bands" . But observe that the first term on the right of Eq. [201 when 
repeated in the whole crystal, diverges because it is Coulomb-like. On the other hand, as 
Fig.[T]shows, the Coulomb tail of the atomic Vg has no importance because the wave function 
never goes far. Then, in using the self-energy potential Vs defined in the atoms we first trim 
the potential with a function as 



Q(r) 
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r > CUT 



r < CUT 
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FIG. 1: (Color online) Self-energy potential (rVs) calculated for the N atom at different ionizations 
/ ranging from 0.5 to —0.2. The lines bunch around that of I = 0.5, which is made thicker. The 
wavefunction for the ionized state (2p) is also shown. Also shown is the potential after the cut-off 
by 9(r). 




The idea is that the atomic self-energy potential is only meaningful where the atomic wave 
function is not negligible. Of course the trimmed self-energy potential is being repeated 
throughout the infinite crystal, so that we are actually calculating "filled hole bands" . Due 
to the trimming, the Coulomb tail (of -1/2 electrons) of the atomic Vs does not penetrate into 
the neighboring atoms. With the trimming, the eigenvalues e a in Eq. [15] become dependent 
on the trimming parameter CUT. However, Eq. [171 sets a recipe to choose the value of CUT: 
one should make the energy gaps extreme. 

The function of Eq. [22] has some important properties: 1 - its derivative is also zero at 
r = CUT, so that its electric field is zero at that point and the cut-off does not add to the 
total charge of the atom; 2 - the trimmed self-energy potential Q(r)Vs(r) is wholly contained 
inside a sphere of radius CUT, which facilitates its use in band-calculation methods such as 
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SIESTA and APW -like. The power n should be large so that the cut-off is sharp. In actual 
practice we tried n = 8 and n = 50 with similarly good results, thus we adopted n = 8, 
which is less abrupt, and does not introduce numerical problems into the programs. Fig. [2] 
shows a typical behavior of a band gap as function of the parameters defining the cut-off 
function. The first increase of the band gap with CUT only means that we are getting 
more of the valence band self-energy as the cut-off is made at larger radii. In principle the 
larger the CUT the more we get of the valence self-energy. After reaching a peak, the gap 
decreases because (i) the potential Vs is penetrating into neighboring atoms, tending to a 
uniform negative potential everywhere in space, shifting downwards all bands, valence and 
conduction alike, and (ii) the self-energy potential perturbation, being broad, diffuses the 
excitation wavefunction thus making it to loose locality and self-energy. In other words, the 
cut-off function should be broad enough so as to include most of the excitation wavefunction 
and thin enough so as not spread it. Thus, the procedure to determine CUT is based in 
Fig. [21 namely we look for the extreme band gap according to Eq. [171 

For a given atom and bonding type, the value of CUT depends little on the chemical 
environment. Because we are using CUT values that make the gaps extreme, small devia- 
tions from the optimal values produce only second-order deviations in the gaps. In the Fig. 
[3] we present the values of anion p state CUT optimized for arsenides, phosphides and ni- 
trides. The anion CUT value has a small dependence on the chemical environment, which 
is approximately linear with the compound bond length. However, the relative variation 
produced in the energy gap is very small. This can be easily verified in the Fig. [2j Around 
the energy gap maximum, the range of CUT found in the optimization for all nitrides leads 
to a change of only 0.05 eV in the energy gap value. This behavior was also verified for all 
other calculated compounds. Therefore, we conclude that is very reasonable to consider the 
same CUT value for the anion potentials. In the Table II I II we show the optimal values of 
the parameter CUT of the trimming function G of Eq. [22j The values for CUT in the Table 
reminds a table of ionic, covalent, or atomic radii, but are not equal to either. 

IV. RESULTS AND DISCUSSION 

We calculated, within the LDA-1/2 approach, the electronic structure for several semi- 
conductors. By comparing the LDA and LDA-1/2 calculation procedures, the LDA-1/2 
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FIG. 2: (Color online) Band gap of GaN as function of the parameter CUT of the cut-off function 
0(r) applied to the self-energy potential of N 2p. The band gap extreme values depend to some 
extent on the exponent n being used (n=8 and 50 are shown). 

calculations lead to no more computational effort than standard LDA, because the values 
of CUT depend only on the atoms, not on their environment, and are calculated just once. 
This is a great advantage of our method. Most of the calculations were made with the 
code VASP "Vienna Ab-initio Simulation Package" using the ultra-soft pseudopotential^ 1 ^ 
to which we added the trimmed self-energy potential. In some instances we repeated the 
calculations with the SIESTA code^ and the results differed by no more than O.leV. The 
two codes are so different, for they use different basis functions and pseudopotentials, that 
the agreement of their results runs in favor of the reliability of our LDA-1/2 procedure. The 
k-space integrals were approximated by sums over a 9x9x9 special-point of the Monkhorst- 
Pack type within the irreducible part of the Brillouin zone.— The number of plane waves for 
the expansion of wave functions was optimized for each system and it is basically the same 
value obtained for optimization of the equivalent standard LDA calculation. The lattice 
parameters used were the experimental ones. 
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The present LDA-1/2 proposal is assuming that in promoting an electron from the valence 
band to the conduction band, the hole thus created is similar to the hole created in the 
atomic photo-ionization. In other words, the hole has the extent of an atomic hole. If the 
hole in the extended system overlapped N equal atoms its self-energy would be 1/N 2 that 
of the atomic hole, and the self-energy potential Vs would be 1/N that of the atom, or it 
would have to be calculated with an ion with 1/2N electrons removed. The results for the 
semiconductors III-V and II- VI, to be presented shortly, definitely point to iV = 1, meaning 
that the hole in the solid resembles much the hole in the atom. In fact, the valence band of 
these semiconductors is known to be made of the anion wavef unctions. On the other hand, 
for the IV elements Ge and Si (and also for diamond), the results point to N = 2 meaning 
that the hole covers the two atoms with covalent bond. 

TABLE III: Values of CUT that make the band gaps extreme, that is, when the self-energy potential 
is defined by Eq. 1181 and trimmed by Eq, [221 The optimal value of CUT, as is the case of an ionic 
or covalent radius, is typical of each atom, and the orbital that was half-ionized. In most cases 
only the anion matters. 
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1.23 
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d 
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The band gaps calculated with LDA-1/2 are presented in Table HVl Here we must remind 
that LDA-1/2 is still a scheme to calculate excitations, not the total energy and the equilib- 
rium lattice parameter. Whereas the LDA results exhibit the well known underestimation 
of the energy gap, LDA-1/2 results present an excellent agreement with experiment. In 
general, by comparing the theoretical LDA and LDA-1/2 band structures we observe that, 
as in LDA-1/2 the self-energy is removed, the valence states are now more localized and are 
pulled down in energy in comparison with the LDA, which results in a larger energy gap. 

The LDA-1/2 entries in Table HVl require a further explanation. In the cases marked 



13 



with two asterisks ** we are adding the trimmed self-energy potential derived from the half- 
ionized anion p-state, and the trimmed self-energy potential derived from the half-ionized 
cation d-state. The questions then are why adding the p-correction to the anion and not 
the s-correction, and why d and not s to the cation. The case of C, Si and Ge, when we 
used -1/4 and not -1/2 ionization, has been discussed above. Thus there seems to be a 
certain degree of arbitrariness in a LDA-1/2 scheme. But that is not so because, from what 
is known from the chemical bonding of these compounds, we could not proceed differently. 
Further we are always keeping in mind the criterion of an extreme band gap (Eq. [T7|) . The 
case of diamond (C) is even more puzzling because we are adding trimmed self-energy s and 
p potentials to a single atom. In this case we are defining the self-energy potential as in 
Eq. [TT]and approaching the method of Filippetti and Spalding^. Again, Eq. [T7] is our guide. 

Figs. HI [5] and [6] depict the corresponding band structures (BS) along the main symmetry 
directions of the Brillouin Zone (BZ) for Si, ZnO, and InN comparing the LDA-1/2 with 
LDA. The zero of energy was placed at the top of the valence band. We chose to show the 
BS for these semiconductors for two reasons: firstly, silicon is one of the most important 
semiconductors, secondly, we would like to show the results for cases where the LDA fails 
completely, as InN, for which LDA gives a semi-metal instead of semiconductor, and ZnO, 
which became a very interesting material with large band gap, and the LDA predicts an 
energy gap much smaller than the experimental value. Moreover, for ZnO it is difficult to 
obtain the correct BS, even performing quasiparticle calculations using GW if the starting 
point is the standard LDA wavef unctions.— From our results, we observe that for Si, the 
LDA-1/2 dispersion relations are similar to the LDA but with the correct band gap energy. 
For the InN and ZnO, the same behavior as Si occurs, but with some differences concerning 
the cation d-states. In both InN and ZnO, the semi core cation d states play an important 
role. In the nitride, the states derived from the atomic 4d (In)-orbital lie close to the 
bottom of the valence band 2s (N)-like and hybridize with fb^l. On the other hand, in 
ZnO, the cation d states lie approximately in the middle of the valence band. Moreover, 
recently it was shown that in both cases the d-states interact and hybridize with the top of 
valence band, and in DFT-LDA there is an underestimation of the binding energies of these 
semicore d states and consequently an overestimation of their hybridization with the anion. 
The enhanced p-d coupling then pushes up the valence-band maximum and the energy gap 
becomes smaller.— i22as&2i By taking these facts into account, for InN and ZnO, LDA-1/2 
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corrects not only the top of valence but also the cation d-states. Thus, by comparing the 
LDA with LDA-1/2 band structures, we observe that the latter has the cation d-orbitals 
deeper in energy This effect is more pronounced in InN. In both cases, ZnO and InN, 
LDA-1/2 is remarkable, leading to values very near the experimental ones. In order to 
study deeply the influence of the cation d-state in both InN and ZnO, in the Table IVl we 
present the results for the valence band width (VBW) and energy gap at different levels 
of the LDA-1/2 procedure. We observe that in the case of ZnO, the O p-state correction 
increases the VBW and the Zn d-state correction decreases the VBW. The combination of 
(0-p) + (Zn-d) corrections presents a smaller VBW than the pure LDA calculation. However, 
both corrections increase the value of the energy gap. The combination of O-p and Zn-d 
corrections results in energy gap in very good agreement with experiment, which again 
states the importance of taking into account the cation d-state. It is worth to point out 
that we obtain this good result, in spite of the fact that the position of Zn d-state (~5 
eV below the top of valence band) is higher in energy than the experimental data (~7.8 
eV).— In the case of InN, as the d state is deeper than in the case of ZnO, the In d-state 
correction is more important for the VBW value, while the N-p correction is more important 
for the energy gap value. Particularly, with (N-p) + (In-d) correction we obtain a value for 
the energy gap which is in good agreement with experiment. Moreover, our full N-p + In-d 
LDA-1/2 calculation is in precise agreement with the measured value obtained from x-ray 
photoemission spectroscopy experiments^, from which the d state of In atom is found to lie 
16.0 eV below the valence band maximum. 

In order to analyze the band dispersion in more detail, we also performed calculations 
to obtain the conduction band effective masses. Thus, now we focus our attention on the 
electronic structure around the conduction-band minima. We fit a parabola to the curves of 
energy versus k around the conduction-band minimum up to 1.0% along the main symmetry 
directions of the Brillouin Zone. Considering the degeneracies and making weighted aver- 
ages we obtain the electron effective masses. Table IVII summarizes the effective conduction 
band masses for several semiconductors. Since a negative value for the LDA-InN band gap 
was obtained, it was not possible to calculate an effective mass in that case and only the 
LDA-1/2 value is shown. We note from the Table that the LDA-1/2 method systematically 
gives larger electron effective masses than LDA. This is due to the fact that with the LDA 
underestimation of band gap energy, the k ■ p interaction between valence (VB) and conduc- 
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tion (CB) bands is stronger leading to smaller effective masses. Therefore, in the cases where 
the correction of the energy gap is more pronounced, the difference between the LDA and 
LDA-1/2 electron effective masses is larger. This is the e.g., the GaAs and ZnO, 

for which the LDA values agree rather poorly with experimental data, and the LDA-1/2 
gives excellent agreement with experiment. Moreover, if we take a look at the whole table 
we observe that the LDA-1/2 effective masses are generally in very good agreement with 
experimental data. 

Therefore, the LDA-1/2 not only improves the band gaps as a "scissor operator" approach, 
but also provides reliable important band structure-derived properties, such as the effective 
masses. 

V. SUMMARY 

The very important problem concerning the calculation of excitations in solids is ad- 
dressed and a method to overcome this problem is developed. The method is inspired in 
the simple half-ionization method. The localization of the hole created by promoting an 
electron from the valence band to the conduction band follows naturally from the method. 
The hole is shown to be representable by a square integrable wavefunction, instead of the 
usual Bloch wave hole of band structure calculations. 

The major success of this method is its reliable description of excited states in solids, 
giving band gap energies, effective masses and band structures in very good agreement with 
experiment, even in the cases for which the LDA markedly fails, e.g., ZnO and InN. The 
method is not more computationally demanding than the LDA calculations. Moreover, the 
method is general and can be applied to a broad class of DFT self-consistent methods, 
all-electron and pseudopotential based. 
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TABLE IV: Band energy gaps (eV) for several semiconductors obtained with the LDA-1/2 at 
experimental lattice constant, by using the VASP code and SIESTA (S), compared with pure 



LDA, GW and experimental results of Ref. |32| . except where noted. Direct energy gaps are noted 
as (d) and indirect as (i). The majority of the LDA-1/2 calculations were obtained using only the 
trimmed self-energy potential of p-anion, exceptions are noted 



a Ref. 



33 





LDA-1/2 


LDA 


Exp. 


GW 


C(i) 


5.25(S)* 


4.13 


5.47° 


5.48-5.77-4 


C(d) 


6.75(S)* 


5.54 


7.3 a 




Si (i) 


1.137, 1.21(S) 0.51 


1.17 a 


1.321&, 0.95-1.104 


Si (d) 


2.9, 2.94(S) 


2.54 


3.05, 3.40 a 




Ge (i) 


0.70 


0.08 


0.66-0.74 a 


0.66-0.834 


A1N (d) 


6.06 


4.27 


6.23 


5.83-6.244 


GaN (d) 


3.52** 


1.95 


3.507 


3.15-3.47 5 


InN (d) 


0.95 ** 


-0.29 


0.7-1.9 


0.20-0.334 


A1P (i) 


2.79 


1.47 


2.52 


2.59 18 


GaP (i) 


2.36(r-L)** 


1.49(r-X) 2.35 


2.55^ 


InP (d) 


1.12** 


0.50 


1.42 


i.Am 


AlAs (i) 


2.73 


1.34 


2.24 


2.15 18 


GaAs (d) 1.41 


0.41 


1.519 


1.231&, 1.40-1.704 


InAs (d) 


0.75 


-0.34 


0.417 


0.3li& 


ZnO (d) 


3.29** 


0.83 


3.4 a 


2.51-3.07-^ 


ZnS (d) 


3.68** 


2.02 


3.91 a 


3.21-3.57-4 



* -l/4p-l/4s 

** -l/2p-anion-l/2d-cation 



Electronic address: guima00@gmail.com 



' Electronic address: mmarques@ita.br 



17 



TABLE V: InN and ZnO valence band width (VBW) and band gap energy values at different levels 
of the LDA-1/2 calculation procedure. The levels presented are: (i) standard LDA calculation, (ii) 
LDA-1/2 anion-p correction only, (hi) LDA-1/2 cation-d correction only, and (iv) full LDA-1/2 
anion-p + cation-d corrections. 

correction VBW (eV) band gap (eV) 



ZnO 



InN 



none 


17.72 


0.83 


O-p 


19.44 


2.14 


Zn-d 


17.01 


1.48 


(0-p)+(Zn-d) 


17.28 


3.29 


none 


15.25 


-0.29 


N-p 


14.49 


1.16 


In-d 


18.23 


-0.49 


(N-p) + (In-d) 


16.85 


0.95 
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FIG. 3: (Color online) The optimized values of the parameter CUT, in which we correct only the 
anion p-state for nitrides, phosphides and arsenides compounds. 
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FIG. 5: Calculated band structures for ZnO (in eV). Dashed lines display LDA and solid lines 
represent LDA-1/2 results. The zero of energy was placed at the top of the valence band. 
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